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Ab initio perturbed ion plus polarization calculations are reported for doubly-charged nonstoi- 
chiometric (MgO)„Mg 2+ (n=l-29) cluster ions. We consider a large number of isomers with full 
relaxations of the geometries, and add the correlation correction to the Hartree-Fock energies for all 
cluster sizes. The polarization contribution is included at a semiempirical level also for all cluster 
sizes. Comparison is made with theoretical results for neutral (MgO)„ clusters and singly-charged 
alkali-halide cluster ions. Our method is also compared to phenomenological pair potential models 
in order to asses their reliability for calculations on small ionic systems. The large coordination- 
dependent polarizabilities of oxide anions favor the formation of surface sites, and thus bulklike 
structures begin to dominate only after n=24. The relative stabilities of the cluster ions against 
evaporation of a MgO molecule show variations that are in excellent agreement with the experimental 
abundance spectra. 
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I. INTRODUCTION 

Cluster science has become a field of intensive research 
both for experimentalists and theoreticians. Small aggre- 
gates have a fundamental interest because they provide 
a link between the molecular and the solid state physics, 
and furthermore they are important in new technological 
applications like nanoelectronics. Elucidating the struc- 
tural and electronic properties of clusters remains a ma- 
jor challenge for present-day science, due to the signifi- 
cant deviations that clusters present in their physical and 
chemical properties when compared both to the molecule 
and the bulk. Another difficulty arises from the huge in- 
crease in the number of different isomers with cluster 
size. Many interesting cluster properties depend largely 
on the cluster structure, at least in the case of ionic and 
covalently-bonded materials, so a complete description of 
the relevant isomer configurations for each cluster size is 
highly desirable. 

In the last years, considerable effort has been de- 
voted to the understanding of metallic and semicon- 
ductor clusters. Meanwhile, studies on metal-oxide 
clusters have been comparatively scarce, in despite of 
their fundamental role in important physical processes 
like heterogeneous catalysis. In this work the inter- 
est is focused on magnesium oxide. Stoichiometric 
MgO clusters have been investigated both experimen- 
tal and theoretically. SaundersoH reported mass spectra 
and collision-induced-fragmentation data for sputter^ 
(MgO)+ cluster ions, and Ziemann and Castlemand' 
performed experimental measurements by rU^ng laser- 
ionization time-of-flight mass spectrornetryBu Theoret- 
ical calculations have been performed at different levels 
of accuracy: simple icmig models based on phenomeno- 
logical pair potpntialsOij were performed by Ziemann 
and CastlemanQ to explain the global trends found in 
their experiments; WilsonEIJ has studied neutral (MgO)„ 
(n<30) clusters by using a compressible-ion modeltil that 
includes coordination-dependent oxide polarizabilities;E3 
semiempirical tight-binding .calculations were reported 
by Moukouri and Noguerajl3o finally, ab initio cal- 
culations on stoichiometric MgO clusters have been 



presented— recently by Recio et al Malliavin and 

CoudrayjI3 and de la Puente et al&3 Nonstoichiometric 
(MgO)„Mg+ and (MgO)„Mg 2+ cluster ions have been 
also detected and, studied experimentally by Ziemann 
and Castleman.Bu££l Specifically, they have obtained, by 
changing the flow rate of the carrier gas in a gas ag- 
gregation source, a mass spectrum comprised almost en- 
tirely of doubly-charged (MgO)„Mg 2+ cluster ions.Ej En- 
hanced stabilities in the small-size regime were found for 
n=8, 11, 13, 16, 19, 22, 25, 27, and were explained in terms 
of compact cubic clusters resembling pieces of the MgO 
crystal lattice. If we exclude the pair potentiaL,calcu- 
lations performed by Ziemann and Castleman,Q'ta there 
is no theoretical investigation of nonstochiometric MgO 
cluster ions. 

The "extra" cation present in nonstoichiometric clus- 
ter ions is expected to result in large structural distor- 
tions whenever a specially compact structure can not be 
constructed. That is not a problem in pair potential cal- 
culations, where the simplicity of the interactions allows 
for a complete geometrical relaxation of the different iso- 
mers at a very modest computational cost. Nevertheless, 
the need to use some selected set of empirical parame- 
ters for Mg 9+ and 9 ~ (q=l,2) ions results in a serious 
questioning of the reliability of these simple ionic mod- 
els. On the opposite side, traditional ab initio calcula- 
tions based on the molecular orbital-linear combination 
of atomic orbitals approximation are very reliable but 
computationally expensive, the computer time require- 
ments scaling with the fourth power of the number of 
atoms in the cluster. Thus, the calculations of Recio et 
al. on (MgO)„ and (MgO)+ clusters£5E3 were performed 
under certain restrictions: (a) the geometries of most of 
the isomers were optimized with respect to a single pa- 
rameter, namely the nearest-neighbor distance; (b) only 
2 or 3 different isomers were considered for each cluster 
size; (c) the correlation energy corrections (calculated at 
the MP2 level) were included onlyJor n < 6. In the calcu- 
lations of Malliavin and CoudrayliZl the geometries were 
more carefully optimized but the results were limited to 
the size range n=l— 6. 

In the present work we report the results of an exten- 
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sive and systematic study of (MgO)„Mg 2+ cluster ions 
with n up to 29. We have employed the ab initio per- 
turbed ion (aiPI) mo del ,E2l which is a Hartree-Fqck model 
based on the theory of electronic separabilityEJEj and 
the_ a6 initio model potential approach of Huzinaga et 
aZ.E3 supplemented with some interaction energy terms 
to account for polarization contributions (see next sec- 
tion). The model has been successfully employed by our 
group injSfiyaral studies of alkali halide clusters and clus- 
ter ions.EfrEj It has been also used in a study, of neu- 
tral stoichiometric (MgO)„ (n=l-13) clusters.tS On one 
hand, our calculations represent a major advance with re- 
spect to pair potential methods, and on the other hand, 
they overcome most of the technical difficulties found in 
more sophisticated ab initio methods: (a) we have al- 
lowed for an appropriate geometrical relaxation of the 
isomers; (b) we have studied a large set of isomers for 
each cluster size (specifically, the total number of iso- 
mers studied is around 400); (c) correlation corrections, 
which have been proved to be essential for an accurate 
description of metal-oxide clustersjia have been included 
for all cluster sizes; (d) we have been able to study rela- 
tively large cluster sizes (up to n=29), enlarging thus the 
usual size range covered by traditional ab initio methods. 

The structural results presented in this work could also 
be useful in the interpretation of possible future exper- 
imental investigations on these clusters. By measuring 
the mobility of cluster ions through an inert buffer gas 
under the influence of a weak electric field, drift tube 
experimental studies ppavide valuable information about 
the cluster geometries Ejtlll 

The rest of the paper is organized as follows: in section 
II we give a brief resume of the aiPI model as applied to 
clusters ( full expositions of the method have been given 
elsewhcrcacj) , a comparison with pair potential models 
which serves to assert the quality of the methodology, and 
the details of the computational procedure. The results 
are presented in section III, and finally section IV sum- 
marizes the main conclusions extracted from this study. 

II. THE AIPI MODEL. COMPARISON TO PAIR 
POTENTIAL MODELS. 

The ab initio perturbed ion mode0 was. originally 
designed for the description of ionic solidsj^j and sub- 
sequentl|V|— adapted to the study of clusters in our 
group.E-3c3 Its theoretical-foundation lies in the theory 
of electronic scparabilityoo and its practical implemen- 
tation in the Hartree-FockJSF) version of the theory 
of electronic separability.EitEa The HF equations of the 
cluster are solved stepwise, by breaking the cluster wave 
function into local group functions (ionic in nature in 
our case). In each iteration, the total energy is mini- 
mized with respect to variations of the electron density 
localized in a given ion. The electron densities of the 
other ions are frozen. In the subsequent iterations each 
frozen ion assumes the role of nonfrozen ion. When the 
self-consistent process (see more details below) finishes, 
the outputs are the total cluster energy E c \ us and a set of 
localized wave functions for each geometrically nonequiv- 



alent ion in the cluster. The cluster .easccy can be written 
as a sum of ionic additive energies o~E3 

N 

Edus = E^ dd , (1) 

where the sum runs over all ions in the cluster, and the 
contribution of each particular ion to the total cluster 
energy (E^ dd ) can be expressed in turn as a sum of in- 
traionic (net) and interionic contributions: 

pi? _ thR , 1 j?RS _ jpR , 1 rpR f \ 

^add — ^net ~T % int ~ net o 

The localized nature of the aiPI procedure has some ad- 
vantages over the usual molecular orbital models. As 
the correlation energy correction in weakly overlapping 
systems is almost intraionic in nature (being therefore 
a sum of contributions from each ion), the localized 
cluster-consistent ionic wave functions may be used to 
attain good estimations of this correction. In this pa- 
per, the correlation energy correction is obt ained through 
dementi's Coulomb-Hartree-Fock method.E3E3 Besides, 
it also allows the development of computationally effi- 
cient codestll which make use of-ihe large multi-zeta ba- 
sis sets of Clementi and RoettO for the description of 
the ions. At this respect, our optimizations have been 
performed using basis sets (5s4p) for Mg 2+ and (5s5p) 
for 2 ~ , respectively. Inclusion of diffuse basis functions 
has been checked and shown unnecessary. Another ad- 
vantage coming from the localized nature of the model 
is the linear scaling of the computational effort with the 
number of atoms in the cluster. This has allowed us to 
study clusters with as many as 59 atoms at a reasonable 
computational cost. 

Selfconsistency has been achieved in the following way: 
for a given distribution of the ions forming the cluster, we 
consider one of them as the active ion R (for instance, a 
particular oxygen anion), and solve the Self-Consistent- 
Field equations for anion R in the field of the remaining 
ions, which are considered frozen at this stage. Next, we 
take another oxygen anion (anion S) as the active ion and 
repeat the same process. If the anion S is geometrically 
inequivalent to anion R, the energy eigenvalues and wave 
functions of electrons in anions S are different from those 
of anions R. We continue this process in the same way 
until all the anions have been exhausted. The same pro- 
cedure is then followed for the magnesium cations. The 
process just described is a PI cycle. We iterate the PI 
cycles until convergence in the total energy of the cluster 
is achieved. Note that the selfconsistent process can be 
accelerated if equivalences between anions or cations are 
imposed by fixing the symmetry of an isomer. In order 
to allow for completely general distortions, we have not 
employed that simplification in the present study. Nev- 
ertheless, equivalences in some ionic wave functions have 
been observed at the end of the calculations for some 
highly symmetrical isomers. 

It is very interesting to compare any quantum model 
designed for the study of ionicjiiaterials with the rigid ion 
and polarizable ion models .lTlj These pair potential mod- 
els are very intuitive, and include in a phenomenological 
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way all the relevant terms in the interatomic potential 
energy. The quality of an ab initio method developed for 
application on ionic materials can thus be assested by 
specifying in which way it improves over a pair poten- 
tial description. The binding energy in a pair potential 
model can be expressed as E bin d = \ v v + J2i E i > 
where Vij is the potential energy between ions i and j, 
and Ef is the self-energy of the ion i measured relative 
to the self-energy of the isolated atom. In the rigid ion 
model, we haved 

yrigid ^m L + Atje-^/P, E?'" sld = 0. (3) 

Tij 

The first term is the electrostatic Coulomb energy be- 
tween two point ions with charges q^ and q_,-, and the 
second term is the short-range repulsive Born-Mayer en- 
ergy reflecting the mutual repulsion due to the overlap 
of the wave functions of the ions. In the polarizable ion 
model, a polarizability an is assigned to each ion so that 
the electron shells can be polarized by the electric field 
created by the other ions in the cluster. Now we haved 

a 2 

ypol _ yrigid . yMD , yDD rpS.pol _ hH_ ,.\ 

V ij V ij " r V V " 1 " v ij ' 2a' 

The self-energy is now different from zero, and the inter- 
action energy contains two new terms: a monopole-dipole 
and a dipole-dipole interaction term. In the Born-Mayer 
repulsive potential, the distance r^ is replaced by an ef- 
fective value to take into account the deformation of 
the electronic shells upon cluster formation. An interme- 
diate pair potential model exists in which the repulsive 
radii of ions are allowed to deform isotropically under the 
effects of other ions in the system, but ionic polarizations 
arising from the positional displacements of each shell 
from its own cope are not considered. This is the breath- 
ing shell modelJ^J and constitutes an important improve- 
ment over the rigid ion model. Improvements over the 
basic polarizable-ion model have also been advanced by 
Madden and coworkersJl3 In the aiPI model, the binding 
energy can be written as a sum of ionic contributions, 
which are in turruaspressed as a sum of deformation and 
interaction termsEfTLJ 

Ebind = ^ Ej% nd = y^(-E£ f + -E^). (5) 

R R 

The interaction energy term is of the form 

T? R — F RS — \^ { pRS , pRS , P RS , P RS \ 

^int — 2-^1 int ~ ^ class +\c + ^ X + ^ overlap ) > 

S^R S^R. 

(6) 

where the different energy contributions are: the clas- 
sical electrostatic interaction energy between point-like 
ions; the correction to this energy due to the finite ex- 
tension of the ionic wave functions; the exchange inter- 
action energy between the electrons of ion R and those 
of the other ions in ther-eluster; and the overlap repul- 
sive energy contribution.Eil The deformation energy term 
Efat is the self-energy of the ion R measured relative to 
the self-energy of the isolated ion. It is an intrinsically 



quantum-mechanical many-body term that accounts for 
the energy change associated to the compression of the 
ionic wave functions upon cluster formation, and incor- 
porates the correlation contribution to the binding en- 
ergy. All those terms are calculated in an ab initio self- 
consistent way. Thus, the improvement over the clas- 
sical rigid-ion description is clear. Improvement over a 
polarizable-ion description can be questionable in princi- 
ple due to one basic assumption used in the actual code, 
namely spherically symmetric electron densities, centered 
on the nuclei. The electron density clouds of the ions are 
thus allowed to distort just isotropically under the ef- 
fect of the other ions in the cluster, and indeed the aiPI 
model could be considered as an ab initio breathing shell 
model. In summary, although the aiPI method treats 
quantum-mechanically all the terms present in a breath- 
ing shell model, plus other many-body terms absent from 
any classical model, it does not describe the dipolar terms 
present in a polarizable-ion model because there are not 
induced pdiuqles. In our previous works on alkali halidc 
clusters, EirEI inclusion of polarization was not considered 
essential. However, the 2 ~ anion has a very deformable 
density cloud (that is a large coordination-dependent po- 
larizability), and dipolar contributions are expected to 
be more important than for halide anions. Relaxing the 
spherical symmetry assumption would allow in principle 
for a proper description of those terms, but many of the 
computational advantages of the aiPI model (which have 
allowed us to perform such a detailed study) would be 
lost. The solution we have chosen is to include the po- 
larization terms in the selfconsistent process with an ex- 
tended polarizable point-ion description.EJ In such a way, 
the "enlarged" model obtained improves clearly over all 
classical descriptions, and can be considered a bench- 
mark for the pair potential calculations. The price to 
be paid is the inclusion of parametrized coordination- 
dependent polarizabilities for the ions. However, the 
polarizabilities that are introduced as parameters in the 
calculations .have been obtained from accurate ab initio 
calculations ,E3 and thus the reliability of the model is re- 
duced just a little compared to the most sophisticated ab 
initio methods, whereas complete geometrical distortions 
can be considered and a large number of isomers studied 
at less cost. Furthermore, the only cases in which doubts 
related to the energetical ordering of the isomers arise 
are those showing near-degenerate isomers, and we will 
show that those situations are not frequent. 

Now we explain the calculational method used in our 
study of (MgO)„Mg 2+ cluster ions. We had no knowl- 
edge a priori as to what shapes can adopt these clusters. 
So we performed an extensive sampling of the potential 
energy surface by generating a large set of random clus- 
ter configurations for each cluster size n with n=l-13. 
Except for the smallest sizes (n=l,2), tipically 100,000 
configurations were generated for each n. All those ran- 
dom configurations were fully optimized by using a rigid- 
ion model. The parameters ia-equation 3 were those used 
by Ziemann and Castleman,cfc3 although some tests wepe 
also performed with the set of parameters of Bush et aljfcH 
Both sets of parameters led to essentially identical struc- 
tures, differing just in minor details. These pair potential 
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calculations could be performed at a low computational 
cost and therefore were very useful in locating reason- 
able initial guesses for the different minima on the po- 
tential energy surface. Then we took tipically the 15-20 
lowest-energy isomers obtained in the pair potential cal- 
culations for each cluster size as input geometries for the 
aiPI calculations. We also studied structures obtained 
by adding or removing in different ways a molecule from 
those pair-potential structures, of course whenever such 
a structure was not already present in the pair-potential 
results. Most of the structures obtained in the range 
n=l-13 could be classified as axbxc, where a, b, and 
c indicate the number of atoms along three perpendic- 
ular edges, so for n>13 we directly studied geometries 
with that formula. We have considered also stackings of 
the different planar structures found for n<13. While 
we can not rule out completely the possibility of hav- 
ing overlooked some important isomer, we believe that 
we have reduced it considerably. We have performed 
aiPI+polarization calculations on all those isomers, with- 
out considering any equivalence between the ions. The 
optimizations of the geometries haj^rbeen performed by 
using a downhill simplex algorithm.Eacj For the oxide pp=. 
larisabilities we have used the values given by Wilson.EiJ 
The magnesium polarizability is taken equal to the bulk 
value .E-J 

We finish this section with a comment about the ionic 
character of our model. Although the aiPI model de- 
scribes the MgO cluster ions in terms of Mg 2+ and 2 ~ 
units, we would like to stress that it does not enter in con- 
flict with any possible assignment of fractional charges 
to each ion in classical models— .In fact, the charges used 
by Ziemann and CastlemarfjEa jsierc +/-1, the charges 
in the potential of Bush et alc3 are +/-2, and both 
potentials give very similar results. The difference in 
the Coulomb part is compensated by a difference in the 
repulsive part. In other words, those charges are just 
parameters. On the other hand, ionic charges can be 
derived from quantuaa-mcchanical calculations following 
the ideas of BaderEJ That has been done recently in 
aiPI calculations on ionic solids-from which fractional 
ionic charges have been derived. EJ Moreover, WilsonEJ 
and Madden and coworkers!!] have shown that an ex- 
tended ionic model including accurate polarization terms 
can account for several effects traditionally atributed to 
"covalency'^ Bulk MgO is excellently described by the 
aiPI modeled, and our aiPI plus polarization results for 
the MgO molecule bond length and vibrational frecuency 
are d= 1.77-Jland ui= 737 cm" 1 , in good-agreement with 
experimenio and ab initio methods l 15 li 17 l The appropri- 
ateness of the method for the study of metal oxide clus- 
ters is thus assested. 



III. RESULTS 

A. Ground State Structures and Low-lying Isomers. 

In the following, we restrict our discussion to the 
structural properties of the lowest-lying isomers. The 
(MgO)„Mg 2+ aiPI+polarization structures are shown in 



Figure 1. The ground state (GS) and one or two low-lying 
isomers are given for each n, except for n>24 where only 
one isomer is shown. A cutoff of 2.54 A, that exceeds by 
« 20 % the Mg-0 distance in the bulk, was arbitrarily 
chosen in order to decide whether an O atom is bonded 
to a Mg neighbor or not. The energy difference with re- 
spect to the most stable isomer is given (in eV) below 
each isomer. For cluster sizes n=l and n=2 (not shown 
in the figure), a chain is obtained as the most stable iso- 
mer. The emergence of the bulk crystalline structure can 
be appreciated from Fig. 1: one-dimensional chains are 
observed for the smallest cluster sizes (n=l-2); a planar 
two-dimensional structure is the GS for n=3; next there 
is a size region (n=4-17) in which three-dimensional (3D) 
structures are already predominant, but without an es- 
tablishment of the bulk symmetry. Specifically, all the 
ions in those structures are in surface-like sites. The 
establishment of face-centered-cubic structures starts at 
n=18, but is not complete until n=24. From that value 
on, all the ground state isomers are fragments of a rock- 
salt lattice. 

Let us describe now in more detail the structures ob- 
tained as a function of n, first for the small cluster sizes 
(n<17). Linear chains are the most stable isomers for 
n=l and n=2. The n=3 GS is a planar rectangular struc- 
ture with a cation attached to a corner. A structure ob- 
tained by removing an anion from a perfect cube appears 
5.12 eV above, and a chain is observed as well as a non- 
competitive isomer. The first three-dimensional ground 
state isomer occurs for n=4. It is obtained by adding 
a cation to a (MgO)4 cubic cluster. A quasiplanar 3x3 
sheet and the chain are not competitive anymore. For 
n=5, a 2x2x3 piece with an anion removed from a cor- 
ner is obtained as the ground state isomer. From n=5 
to n=17, all the ground state isomers are obtained by 
adding a MgO molecule either to the GS or to a low- 
lying isomer of the (MgO) n _iMg 2+ cluster, in such a 
way that no bulk ions (ions with coordination six) are 
present. The formation of surface sites (with coordina- 
tion equal to 3, 4, or 5) is favored instead. This is a 
direct consequence of the larger polarizabilities of oxide 
anions with lower coordinations. The minimization of 
just the electrostatic Coulomb energy between point ions 
would lead to the formation of more compact bulklike 
structures, for which the ions would tend to attain its 
full first-coordination sphere. But when polarization ef- 
fects are included, the dipole stabilization energy has to 
be minimized also, favoring a reduction of the coordina- 
tion number and the formation-af surface sites, the same 
conclusion achieved by Wilsor£3 in his study of neutral 
(MgO)„ clusters. As a result that is really worth of men- 
tioning, the GS structure for n=13 is not a 3x3x3 perfect 
cube, that is what one would have "expected" by compar- 
ing with t he situation encountered for alkali halide cluster 
ionsOH'Lil This essentially different structural behavior 
has its roots in the increased polarizability of the oxide 
compared to halide anions, and can be accounted for by 
an enlarged ionic model like the one used here., without 
resorting to the inclusion of "covalent" effects. EJ'Ej The 
3x3x3 (MgO)i3Mg 2+ isomer has an oxide anion with 
coordination 6, so with a small polarizability. Moreover, 
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the corner sites of that isomer are occupied by cations, 
so there is not any oxide anion with coordination 3. All 
this results in a small polarization contribution to the en- 
ergy, that is not compensated by the increased Madelung 
energy. 

Eventually there has to be a cluster size where bulklike 
fragments begin to dominate. The situation changes a lit- 
tle at n=18, where a ground state isomer with bulklike 
ions appears for the first time, specifically a 3x3x4 cubic 
structure with a cation attached to a corner. Let us note 
that a major difference between this structure and the 
3x3x3 (MgO)i3Mg 2+ isomer is that now there are oxide 
anions in corner positions which enhance the polariza- 
tion contribution. This enhancement, together with the 
increased Madelung term associated to a bulklike frag- 
ment, make that isomer more stable than the one ob- 
tained by adding a MgO molecule to the (MgO)iyMg 2+ 
ground state. The establishment of bulklike symmetry is 
not complete yet, however. For n=19, a 3x3x4 struc- 
ture with three ions attached is still favored, even more 
than for n=18, because now none of the four oxide an- 
ions with coordination 3 is capped by a magnesium. By 
adding 5 ions to such a cubic structure, however, at least 
two oxide anions in corner sites are capped, and so the 
GS structure of (MgO)2oMg 2+ is again an isomer with 
only surface sites. For n=22, a 3x3x5 perfect structure 
can be formed, in analogy with the case of n=13. Such a 
structure would have no oxide anions in corner positions 
and two oxides with coordination 6, however, and it is not 
still energetically favored. From n=24 on, all the struc- 
tures obtained are fragments of a bulklike crystalline lat- 
tice. Although it is possible that for some larger cluster 
size the ground state structure will be not yet the most 
compact face-centered-cubic fragment, we think we can 
accept with relative safety n=24 as the critical size where 
bulklike symmetry emerges. This is consistent with the 
critical size (n=30) estimated by WilsonE] for neutral 
(MgO) n clusters. 

Our conclusions suggest substantial structural differ- 
ences between (MgO)„Mg 2+ cations and (MgO)„0 2 ~ an- 
ions. For example, the perfect 3x3x3 cube isomer of 
(MgO)i30 2 ~ would not have any oxide anion with co- 
ordination six, and all the eight corner sites would be 
occupied by oxides. The case of 3x3x5 (MgO)22 2 ~ 
isomer would be similar. We have performed additional 
calculations for those two cluster sizes and have obtained 
the 3x3x3 and 3x3x5 fragments as ground state isomers 
of (MgO)i 3 2 - and (MgO) 22 2 -, respectively. 

A comparison with the results obtained for .alkali halide 
cluster ionscZrEa and neutral (MgO) n clustersEJ'ta is illus- 
trative. Hexagonal prysm atip structures are very com- 
mon for (MgO)„ clusters, EEnl3-as-4vell as for some neutral 
alkali halide (AX) n clusters. Ejl3 However, cubic struc- 
tures are clearly predominant for (MgO)„Mg 2+ cluster 
ions, being this the main effect of nonstoichiometry and 
net charge: (a) perfect cubic structures can be built up 
with an odd number of ions, but an even number of ions 
is needed to construct a perfect hexagonal prysm; (b) for 
"defect-like" structures, the extra charge is screened more 
effectively in the cuboid-like isomers. Examples pan be 
found at n=3, where the hexagonal (MgO) 3 isomerE3 con- 



verts into a rectangular structure upon adding Mg 2+ to 
lower the repulsion between cations, and at n=10, where 
adding a triatomic Mg-O-Mg to the hexagonal prysm of 
the neutral (MgO^Ej results in an isomer less stable than 
that obtained by adding Mg-O-Mg to a rocksalt piece 
(second isomer in figure 1). Let us now compare wthi 
the results obtained for alkali halide cluster ions£3€3 
The main difference comes from the polarization contri- 
bution. As stated above, that is less important in alkali 
halide materials, for which the establishment of bulklike 
symmetry is already complete at n=13. There is just an- 
other minor difference: centered structures, containing 
an inner alkali cation with high coordination, were found 
frecuently for (NaI) n Na + and (CsI)„Cs + , but such iso- 
mers are not observed for (MgO)„Mg 2+ . Evidently, the 
repulsion between oxide anions is larger than the repul- 
sion between halide anions, and allocating six 2 ~ anions 
around a Mg 2+ cation in energetically less favorable. 

We finish this section with a discussion of the validity 
of pair potential models and the effects of the polariza- 
tion correction on cluster structure. Comparison of our 
pair potential calculations with the aiPI results lead us 
to the following two important points: (a) pair potential 
models fail to reproduce quantitatively the expansion of 
the interionic distances with cluster size; (b) the energet- 
ical ordering of the different isomers for a given cluster 
size is not reliable. On the positive side, the shape of the 
isomers is usually quite well reproduced, and the main 
structural correction obtained with an aiPI calculation 
is a global scaling of all the interionic distances. This 
last observation supports the use of pair potentials to lo- 
cate good initial configurations for ionic clusters, from 
which ab initio calculations can be started. The effects 
of the polarization correction are sizable for n< 17. The 
correct determination of the ground states would not be 
achieved with a model that did not include polarization 
corrections, as surface sites would not be specially fa- 
vored. Moreover, for the smallest cluster sizes (n<8), the 
dimensionality of the isomers is not reproduced either. 
Specifically, linear chains are the GS isomers for n=l 
3 and planar structures are the most stable isomers for 
n=4-7 in aiPI calculations without polarization. Thus, 
inclusion of polarization speeds up the emergence of 3D 
structures. The imprecisions associated to the lack of 
polarization corrections turn smaller as the cluster size 
increases and bulklike fragments begin to dominate. 

B. Relative Stabilities. Comparison to Experimental 
Results. 

In order to study 

the relative stabilities of (MgO)„Mg 2+ cluster ions we 
calculate the evaporation threshold energydO required 
to remove a MgO molecule from the ground state iso- 
mer. For a cluster (MgO)„Mg 2+ this is done as follows: 
we consider the optimized GS structure and identify the 
MgO molecule that contributes the least to the cluster 
binding energy. Then we remove that molecule and relax 
the resulting (MgO)„_iMg 2+ fragment to the nearest lo- 
cal minimum. The total energy required to evaporate a 
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molecule is then the difference between the energy of the 
parent cluster and the energies of the two fragments, one 
of which is the MgO molecule: 

Eevaporationi^) — ^cluster 

+E(MgO) - E cluster [(MgO) n Mg 2 +}. (7) 

This process can be termed locally adiabatic because 
both fragments are allowed to relax to the local minimum 
energy configuration after the evaporation. For some 
cluster sizes, the fragment of size (n-1) left when a MgO 
molecule is removed from (MgO)„Mg 2+ does not lie on 
the catchment basin of the (MgO)„_iMg 2+ GS isomer, 
so that the evaporation threshold energies are larger than 
the energy differences between adjacent ground states mi- 
nus E(MgO) in those cases. The evaporation energies 
are plotted as a function of n in figure 2. In the ex- 
periments performed by Ziemann and Castlemanjla the 
abundances in the mass spectra reflect the relative sta- 
bilities of (MgO)„Mg 2+ cluster ions against evaporation 
of a MgO molecule. Starting from the largest cluster 
size studied, we can appreciate minima in figure 2 for 
n=28, 26, 23, 20, 17, 14, 12, 9, and 5. Thus, evaporation 
of a MgO molecule from those clusters is easy relative 
to evaporation from clusters of other sizes. Moreover, 
the curve shows abrupt increases precisely at those clus- 
ter sizes, so that evaporation of a MgO molecule from 
the clusters of size (n-1) is significantly more expensive 
than evaporation from cluster ions of size n. All of this 
results in an enrichment of (MgO)„Mg 2+ clusters with 
n=4, 8, 11, 13, 16, 19, 22, 25, and 27, a result which is in 
complete agreement with the experimento Nevertheless, 
there is not a total correspondence between the experi- 
mental magic numbers and maxima in the evaporation 
energy curve. The experimental enhanced abundances 
result from a balance between two main processes: (a) 
evaporation from clusters with size (n+1) results in an 
enrichment of the (MgO)„Mg 2+ clusters; (b) evaporation 
from clusters with size n has the opposite effect. A quan- 
tity that measures the stability of a (MgO)„Mg 2+ cluster 
ion relative both to (MgO)„+iMg 2+ and (MgO)„_iMg 2+ 
cluster ions is thus 

^2(^7 — Ej eva p 0ra ij i0n {jl -\- 1) -E eva p 0ra fa on (Tl') , (8) 

where both terms in the difference are calculated using 
eq. (7) as explained above. This quantity is plotted as a 
function of n in Figure 3. Minima are apparent at n=4, 
8, 11, 13, 16, 19, 22, 25, and 27. Thus, the experimental 
results concerning enhanced relative stabilities are prop- 
erly reproduced by our calculations. ■— ■ 

Opposite to the case of alkali halide clustersc3 and clus- 
ter ions,C3 the magic numbers observed for (MgO) n Mg 2+ 
can not be clearly explained in terms of a compactness 
argument. Instead, all the magic numbers observed ver- 
ify the following: when removing the least bound MgO 
molecule from (MgO)„Mg 2+ (where n is a magic num- 
ber), the resulting (MgO)„-iMg 2+ is never in the catch- 
ment area of the ground state isomer, and thus the evap- 
oration energy is large. Let us note that in some cases 
(n=13,19,25) this is not a direct consequence of the shape 
of the ground state isomers involved, but of the asymme- 
try between the different facets of those isomers, which 



make energetically more favorable to allocate three "ex- 
tra" ions on a different facet than one "extra" ion. More- 
over, the magic numbers (n=4, 11, 13, 16, 22, 25, 27) verify a 
second property: when removing the least bound MgO 
molecule from (MgO)„ + iMg 2+ (where n is again a magic 
number), the resulting (MgO)„Mg 2+ structure is either 
the ground state structure or a low-lying isomer, and thus 
the evaporation energy shows a local minimum for those 
cluster sizes. Although there is a striking periodicity of 
three in the magic numbers observed between n=13 and 
n=27, our results indicate that there is not any special 
structural motif behind that periodicity. 

IV. SUMMARY 

The ab initio perturbed ion model, supplemented with 
a semiempirical treatment of dipolar terms, has been em- 
ployed in order to study the structural and energetic 
properties of (MgO)„Mg 2+ (n=l-29) cluster ions. We 
have used a rigid ion model to generate a large set of 
initial configurations for each cluster size. Next, those 
structures have been optimized without consider any 
equivalence between different ions, allowing thus for an 
appropriate geometrical relaxation. Correlation correc- 
tions, which turned out to be-essential for a proper de- 
scription of (MgO)„ clustersjla have been added to the 
Hartree-Fock energies for all cluster sizes. Inclusion of 
parametrized coordination-dependent polarizabilities re- 
duces just a little the reliability of the model compared 
to full ab initio methods, because those polarizabilities 
have been themselves extracted from accurate ab ini- 
tio calculations, and allows to perform such a system- 
atic study with appropriate geometrical relaxations at a 
significantly reduced computational cost. The structural 
trends of (MgO)„Mg 2 + cluster ions have been described. 
The emergence of the bulk crystalline symmetry proceeds 
in different stages, from one-dimensional configurations 
for n=l-2 to a planar structure for n=3, and then to 
three-dimensional structures with a clear predominance 
of cubic structures with surface sites for n=4-17. The 
critical size at which bulklike rocksalt fragments dom- 
inate the spectrum is n=24. For (MgO)i3Mg 2+ and 
(MgO) 22 Mg 2 +, the "expected" 3x3x3 and 3x3x5 rock- 
salt isomers are neither the ground states nor low-lying 
isomers. This at first sight surprising result has been 
rationalized in terms of the importance of the polariza- 
tion contribution, and the preference of oxide anions for 
surfacelike (specially corner) sites. A comparison .with, 
the structures obtained for neutral (MgO)„ clusteralj'Ea 
shows that the preference for cuboid-like structures is a 
direct consequence of nonstoichiometry and net charge. 
Comparison_with the structures obtained for alkali halide 
cluster ionsEfl reveals substantial differences that can be 
explained in terms of the larger polarizabilities of oxide 
anions compared to halide anions, and to a less extent in 
terms of the larger repulsion between oxide anions com- 
pared to the repulsion between halide anions. The main 
effects of the polarization correction is to speed up the 
emergence of 3D structures, to lower the average coordi- 
nation, and to favor the formation of surface sites. Once 
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bulklike structures are well stablished polarization plays 
a less important role. 

The relative stabilities of (MgO)„Mg 2+ cluster ions 
have been studied by calculating the energy required to 
evaporate a MgO molecule from the GS structures and 
the first differences between those evaporation energies, 
that reflect the stability of a cluster of size n relative 
to the stabilities of the neighbor clusters of sizes (n+1) 
and (n-1). Both sets of calculations predict the clus- 
ters with n=4, 8, 11, 13, 16, 19, 22, 25, and 27 to be 
specially stable, a result that is in complete agreement 
with the experimental enhanced abundances reported by 
Ziemann and Castleman.Eia The interpretation of the en- 
hanced stabilities in terms of highly compact fragments 
of a face-centered-cubic crystalline lattice is not appro- 
priate for such small clusters. Those enhanced stabilities 
have been explained in terms of an analysis of the evap- 
oration process, which involves the explicit consideration 
of isomer structures different from the ground states. 

Finally, a comparison has been made with the re- 
sults obtained by using a phenomenological pair poten- 
tial model. This comparison has shown that, while pair 
potential calculations are really helpful when obtaining 
initial guesses for the isomer geometries, the distances 
and the energetical ordering of the isomers obtained with 
them are generally not reliable. 
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Captions of figures 

Figure 1. Lowest-energy structure and low- lying iso- 
mers of (MgO)„Mg 2+ cluster ions. Light balls are Mg 2+ 
cations and dark balls are O 2- anions. The energy dif- 
ference (in eV) with respect to the most stable structure 
is given below the corresponding isomers. For n>24 only 
the ground state isomer is shown. 

Figure 2. Evaporation threshold energies required to 
remove a neutral MgO molecule from (MgO)„Mg 2+ clus- 
ter ions as a function of n. The local minima in the 
evaporation energy curve are shown explicitely. 

Figure 3. Second derivative A 2 (n) for (MgO)„Mg 2+ 
cluster ions as a function of the cluster size. Minima 
identify those cluster sizes with an enhanced relative sta- 
bility. 
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